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Abstract-  The  additional  surface  roughness  created  by 
drops  impacting  the  sea  surface  can  bias  wind  speed  estimates 
obtained  from  satellite  scatterometric  measurements.  The 
additional  roughness  essentially  depends  on  the  rain  content 
in  very  large  drops.  The  estimate  of  this  contribution  will  be 
highly  dependent  on  the  model  chosen  for  the  drop  size 
distribution.  However,  it  has  been  observed  that  the 
reflectivity  of  the  drops  falling  in  the  atmosphere  has  a  drop- 
size  dependence  that  is  very  similar  to  that  of  the  surface 
roughness  induced  by  rain.  This  paper  shows  how  reflectivity 
data  can  be  used  to  improve  the  estimate  of  the  elevation 
variance  resulting  from  drop  impact. 

I.  INTRODUCTION 

The  impact  of  raindrops  on  the  sea  surface  can  significantly 
modify  its  reflectivity.  This  can  introduce  a  bias  in  the 
wind  speed  estimations  based  on  scatterometric  or 
radiometric  measurements.  Moreover,  as  explained  by 
Meneghini  et  al  [9,10],  rain  rate  retrieval  algorithms 
generally  also  rely  on  the  estimation  of  the  surface 
reflectivity.  The  latter  is  often  considered  equal  to  the 
reflectivity  just  outside  the  rain  cell,  which  does  not 
account  for  the  effects  of  drop  impact  and  could  introduce  a 
corruption  of  the  estimated  rain  attenuation  appearing  in 
retrieval  algorithms  [5]. 

As  a  first  step,  we  assume  that  the  impact  of  drops  on  the 
sea  surface  yields  an  additional  component  in  the  surface 
elevation  spectrum,  superimposed  to  the  surface  spectrum 
due  to  instantaneous  wind  speed  and  wind  history.  In  order 
to  quantify  this  effect,  laboratory  experiments  were  carried 
out  by  Bliven  and  Sobieski  [1,2,11,12].  Those  yield  the 
roughness  appearing  on  a  water  surface  subject  to  artificial 
rain,  along  with  corresponding  radar  backscattering 
signatures.  The  two  spectra  are  represented  in  Fig.  1  for 
increasing  wind  friction  velocities  (u*)  and  rain  rates  R. 
The  wind  spectrum  is  obtained  here  from  the  model  of 
Bjerkaas  and  Riedel  (1979)  but  could  be  replaced  by  any 
other  one  [8]. 

Two-scale  electromagnetic  scattering  models  were  then 
used  to  quantify  the  effect  of  rain  on  the  surface  reflectivity 
change.  In  some  configurations,  it  has  been  found  that  the 
surface  radar  reflectivity  due  to  rain  could  change  by  up  to 
10  dB  [3]. 


II.  SURFACE  ELEVATION  SPECTRUM 

The  wave  tank  experiments,  associated  with  scattering 
simulations,  were  performed  initially  for  mono-disperse 
artificial  rain  events,  only  made  of  drops  with  diameters 
equal  to  2.8  mm.  However,  natural  rain  being  far  from 
mono-disperse,  further  experiments  to  analyze  the  drop- 
size  dependence  were  conducted  [7]. 
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Fig.l:  Radial  sea  surface  elevation  spectrum  S(k)  for 
increasing  wind  friction  velocities  u*  and  rain  rate  R.  u* 
ranges  from  20  to  120  cm/s  by  steps  of  10  cm/s  and  R 
ranges  from  20  to  100  mm/hr  by  steps  of  20  mm/hr. 

From  those  new  experiments,  it  appeared  that  the  largest 
drops  have  a  much  larger  contribution  to  the  surface 
elevation  variance  than  small  ones.  This  was  consistent 
with  single-drop  experiments  performed  by  Craeye  et  al. 
[4],  who  found  that  the  energy  transferred  by  the  drop  into 
surface  waves  is  not  simply  proportional  to  the  kinetic 
energy  of  the  drop. 

Indeed,  the  larger  the  drop,  the  larger  the  relative 
contribution  to  the  surface  energy.  Actually,  it  appeared 
that  the  contribution  of  a  given  drop  to  the  surface 
elevation  variance  is  practically  proportional  to  its  squared 

2  2 

momentum  m  V  .  This  simple  rule  seems  to  be  consistent 
with  the  experimental  data  analyzed  for  various  mono- 


0-933957-31-9 


528 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

01  SEP  2003 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

A  relationship  between  atmospheric  rain  reflectivity  and  elevation 

5b.  GRANT  NUMBER 

viuimitc  uuc  iu  uiup  uupiiti  uii  me  suimtc 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Universite  catholique  de  Louvain  Place  du  Levant,  2,  B- 1348 
Louvain-la-Neuve,  Belgium 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM002146.  Oceans  2003  MTS/IEEE  Conference.  Held  in  San  Diego,  California  on  September 
22-26,  2003.  U.S.  Government  or  Federal  Purpose  Rights  License,  The  original  document  contains  color 
images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

4 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


disperse  rain  events,  as  reported  in  [7].  Figs.  3  and  4 
illustrate  the  quality  of  the  scaling  rule  based  on  the  square- 
momentum  law. 


Fig.  2:  Surface  elevation  variance  obtained  in  the 
laboratory  for  different  drop  diameters. 


Fig. 3.  Surface  elevation  variance  rescaled  by  the  flux  of 
squared  drop  momentum.  From  the  2.3  and  2.8  mm  cases 
toward  the  (unchanged)  4.2mm  case. 

Knowing  that  only  the  largest  drops  will  contribute  to  the 
surface  energy  associated  with  rain,  we  will,  in  a  first  stage, 
omit  the  dependence  of  fall  velocity  on  the  drop  size.  This 
simplification  is  justified  by  the  fact  that,  for  drops  larger 
than  3  mm,  the  error  introduced  on  the  velocity  is  of  a  few 
percent  only  [6].  The  dependence  of  velocity  on  drop  size 
will  be  reintroduced  in  the  next  sections. 

Hence,  concentrating  on  the  mass  appearing  in  the  square 
momentum  law,  we  observe  that  the  contribution  of  a  drop 
to  the  surface  energy  is  proportional  to  the  sixth  power  of 
the  drop  diameter.  Consequently,  in  the  case  of  natural  rain, 
the  estimate  of  the  surface  elevation  variance  will  strongly 
depend  on  the  choice  made  for  the  drop  size  distribution 
(DSD). 


In  the  following,  we  show  that  this  difficulty  can  be 
overcome  by  directly  linking  the  surface  elevation  variance 
with  the  reflectivity  of  the  drops  in  the  atmosphere.  The 
reason  is  that,  in  the  low-frequency  approximation,  the 
contribution  of  a  drop  is  also  proportional  to  the  sixth 
power  of  the  drop  diameter  [14].  In  other  terms,  the 
reflectivity  of  the  drops  can  be  regarded  as  a  good  image  of 
the  perturbation  they  cause  on  the  sea  surface,  and  the 
relationship  between  these  two  quantities  should  be  almost 
independent  from  the  model  chosen  for  the  drop  size 
distribution.  In  the  following,  we  will  illustrate  this 
relationship  by  taking  into  account  18  different  DSD's.  In 
this  exercise,  we  will  assume  the  drops  to  be  spherical,  and 
the  reflectivity  data  are  produced  from  exact  Mie 
computations,  assuming  lossy  dielectric  spheres.  In  the 
following,  the  models  for  the  surface  elevation  variance 
and  drop  reflectivity  will  be  recalled,  and  they  will  be 
confronted  with  each  other. 


II.  SURFACE  ELEVATION  VARIANCE  AND 
RAINDROP  REFLECTIVITY 


Fig.4.  Normalized  radar  cross-section  of  a  spherical  water 
drop  at  13.8  GHz. 

The  energy  contributed  by  a  drop  to  the  surface  lasts  for  a 
certain  time,  which  we  will  call  the  "lifetime"  of  the  ring- 

waves  Tx .  For  natural  rain,  the  surface  elevation  variance 

2 

G  is  obtained  by  averaging  the  ring-wave  energy  over 
the  drop  size  distribution  n(D).  It  yields: 

cr=T|  (/?)  ja,  D6v2(D)v(D)n(D)  dD  (1) 

0 

where  ax  is  a  constant,  and  the  lifetime  Tx  is  assumed  to 
be  depend  on  the  rain  rate  R  only.  Its  value  depends  on  the 
interactions  of  ring-waves  with  each  other,  and  with  the 
other  events  resulting  from  drop  impact  (splashes, 
turbulence).  For  light  rain  rates,  the  lifetime  is  constant, 
because  the  interactions  are  negligible,  while  it  tends  to 
become  smaller  when  the  rain  rate  increases.  tx(R)  is 
obtained  within  a  multiplicative  constant  from 
experimental  data  obtained  with  mono-disperse  rain  [3], [7]. 

The  reflectivity  from  the  drops  in  the  atmosphere  can  be 
obtained  using  the  Mie  scattering  theory.  For  drops  that  are 
very  small  compared  to  the  wavelength  (at  least  ten  times 
smaller),  the  backscattering  cross-section  is  proportional  to 
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the  sixth  power  of  the  drop  diameter  [10].  As  explained 
above,  this  explains  the  functional  similarity  between  the 
surface  elevation  variance  and  the  atmospheric  drop 
reflectivity. 

z  =  f  f(D)n(D)dD  with  /(D)  *  D6  (2) 

Jo 

However,  the  condition  referred  to  above,  which  stands  for 
the  Rayleigh  scattering  approximation,  is  not  necessarily 
fulfilled  at  the  frequencies  used  in  existing  satellite 
configurations.  Hence,  exact  Mie  computations  [14]  have 
been  carried  out  here,  assuming  spherical  drops.  Fig.  2 
shows  the  radar  cross-section  of  a  drop  obtained  versus 
drop  diameter  at  13.8  GHz,  which  is  the  frequency  of  the 
TRMM  precipitation  radar  [10].  At  that  frequency,  the  drop 
complex  dielectric  constant  is  close  to  36  -  35  j  [13]. 
Obviously,  the  radar  is  sensitive  to  the  very  largest  drops 
only. 


Fig. 5.  Surface  elevation  variance  computed  with  the  help  of 
18  different  drop  size  distributions  found  in  [10]. 


III.  EFFECTS  OF  THE  VARIABILITY  OF  DSD 
MODELS 

The  drop  size  distribution  (DSD)  models  reported  in  the 
literature  are  widely  deviating.  In  this  paper,  we  selected 
18  of  the  22  distributions  gathered  by  Montanari  [11].  Four 
distributions  have  been  discarded  because  of  their  poor 
rain-rate  consistency.  Using  the  notations  defined  above, 
the  rain  rate  should  be  equal  to: 

oo 

R  =  0.6  10  3  j n(D,  R)D3v(D)dD  (3) 

0 

Figs.  5  and  6  show  the  surface  elevation  variance  and 
reflectivity  factors  obtained  with  the  help  of  these  DSD's. 
The  uncertainty  resulting  from  the  choice  of  the  DSD 
model  is  of  the  order  of  a  6:1  ratio.  As  it  is  not  possible  to 
say  a  priori  which  is  the  best  DSD  model  (the  actual  DSD 
depends  on  many  factors,  like,  for  instance,  the  type  of 


clouds),  it  is  very  difficult  to  obtain  an  estimate  of  the 
surface  energy  induced  by  impacting  drops.  It  should  be 
noted  that  some  DSD's  lead  to  very  high  values  of  the 
elevation  variance.  This  is  due  to  the  fact  that  many  DSD's 
over-estimate  the  large-drop  density. 


Fig. 6.  Reflectivity  factor  computed  with  the  help  of  18 
different  drop  size  distributions  found  in  [11]. 

It  is  interesting  to  realize  that  the  integral  I  appearing  in  the 
expression  of  the  elevation  variance  (1)  is  similar  to  the  one 
appearing  in  the  reflectivity  factor  (2).  Hence,  those  two 
quantities  are  almost  proportional  to  each  other.  In  other 
terms,  the  relationship  between  drop  reflectivity  and 
elevation  variance  presents  a  weak  dependence  on  the 
actual  model  chosen  for  the  DSD.  This  is  visible  in  Fig.  7, 

where  the  dashed  lines  correspond  to  the  use  of  a  D6 


Fig  7.  With  the  help  of  the  same  DSD  models  as  before,  I 
integral  appearing  in  the  expression  of  the  surface  elevation 
variance,  versus  the  reflectivity  factor.  The  dashed  lines  are 
obtained  when  the  the  Rayleigh  approximation  is  used. 
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model  for  the  drops  cross-section,  instead  of  the  exact  Mie 
computations.  Obviously,  when  the  surface  elevation 
variance  is  obtained  from  reflectivity  data,  the  choice  of  the 
DSD  model  becomes  a  much  less  critical  question. 

Fig.  7  only  shows  the  value  of  the  integral  appearing  in  the 
expression  of  the  surface  elevation  variance.  For 
completeness,  the  lifetime  function  should  also  be 
estimated.  We  assumed  that  it  depends  only  on  the  rain 
rate.  In  most  cases,  there  is  no  independent  data  available 
on  the  rain  rate,  such  that  it  has  to  be  delineated  from  the 
drop  reflectivity.  As  the  rain  rate  (3)  and  the  drop 
reflectivity  (2)  do  not  have  similar  drop- size  dependences, 
the  estimated  rain  rate  will  be  strongly  dependent  on  the 
choice  of  the  DSD.  This  constitutes  the  standard  challenge 
of  rain  rate  estimation.  However,  it  should  be  recalled  that 
the  lifetime  x  has  a  weak  dependence  on  the  rain  rate. 
Hence,  the  presence  of  the  lifetime  in  the  expression  of  the 
surface  elevation  does  not  dramatically  enhance  the 
dependence  of  the  result  on  the  DSD.  This  is  shown  in  Fig. 
8. 


Fig. 8.  With  the  help  of  the  same  DSD  models  as  before, 
surface  elevation  variance  versus  reflectivity  factor. 

IV.  CONCLUSION 

Rain  can  strongly  affect  the  reflectivity  from  the  surface  of 
the  sea.  The  effect  of  the  impacting  drops  on  the  surface 
roughness  increases  very  fast  with  drop  size.  As  a  result, 
only  the  largest  drops  affect  the  surface  roughness.  This 
effect  is  very  difficult  to  estimate  when  only  the  rain  rate  is 
known,  because  the  large  drop  content  strongly  depends  on 
the  models  chosen  for  the  DSD  (cf.  Fig.  5).  This 
uncertainty  can  be  dramatically  reduced  when  drop 
reflectivity  data  are  available,  as  it  is  the  case  for  the 
TRMM  data.  Indeed,  drop  reflectivity  and  surface  elevation 
variance  correspond  to  almost  identical  moments  of  the 
DSD.  This  is  illustrated  in  Fig.  8,  where  we  see  that  the 
about  6:1  uncertainty  is  now  reduced  to  about  2:1.  The 
remaining  uncertainty  is  caused  by  (i)  the  departure  of  the 
Mie  solution  from  the  Rayleigh  approximation  at  high 


frequencies  and  (ii),  the  dependence  on  the  surface 

elevation  on  the  lifetimes  of  the  ring- waves.  The  latter 

quantity  certainly  deserves  further  attention. 
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